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■ Two-neutron correlation in the ^ S channel in quasi two-dimensional (2D) neutron systems at zero 
(N . temperature is studied by means of the BCS theory with finite-range effective nuclear forces. The 

$H ' dineutron correlation in low density neutron systems confined in an infinite slab is investigated in 

a simplified model that neutron motion of one direction is frozen. When the slab is thin enough, 
two neutrons form a tightly bound dineutron with a small size in the quasi-2D system, and a Bose 
dineutron gas is found in low density limit. With increase of Fermi momentum, the neutron system 
^ : changes from the Bose-gas phase to the superfiuid Cooper pair phase. The density dependence of 

the neutron pairing shows the BCS-BEC crossover phenomena at finite low-density region. In the 
transition region, the size shrinking of neutron pair and enhancement of pairing gap are found. The 
relation to dineutron correlation at surface of neutron-rich nuclei is also discussed. 

'b 

I K I. INTRODUCTION 

(N ' 

^ ■ Two-neutron correlations in neutron-rich nuclei is presently one of the fore-front subjects in the physics of unstable 
[ — ] nuclei. In two-neutron halo nuclei [l[ such as ^^Li, the dineutron correlation was theoretically predicted in many 
T-H studies 0, H, 0, [il, |7| |8j, |9|, [l^, [ll| and is supported by experiments [H, [H, [13, [H, [11] ■ Dineutron correlations are 

■ discussed also in light neutron-rich nuclei such as in ®He [ij, [Tsl. [l9l a nd also in medium-heavy neutron-rich nuclei 
, [HjUlj as well as asymmetric nuclear matter (for example, Refs. |22l. |23|. [2^. [25l. [26l. [27l. [28l [29t and references therein). 

' Dineutron correlations are characterized by strong spatial correlations of two neutrons in the channel. Although 
, the neutron- neutron interaction is too weak to form a two-body bound state, the attraction in the channel is 

■ rather strong as is known from the large scattering length as = — 18.5± 0.4 fm [s^l- Originally the possible existence 
of a dineutron near the surface of nuclei was predicted by Migdal in 1972 [3l|. In his work, it was shown that two 
interacting particles in a potential well can form a bound state if there is a single-particle level with energy close to 
zero. This idea describes the binding mechanism of a Borromean system and is extended to the dineutron cluster 

^ ^ , picture and the dineutron correlation in two-neutron halo nuclei. 

Recently, dineutron correlations have also been discussed from the point of view of Bosc-Einstein condensation 
(BEG), which is considered along similar lines as deuteron and a-particle condensation suggested in dilute nuclear 
matter [2^ [H, [H, [H, HH] as well as a condensate states in excited states oi Z = N = even nuclei [H, 113, HI] . This 
nn-BEC feature stems from the surprisingly small size of Cooper pairs recently revealed in the surface of nuclei [2l[ . 
In nuclear systems, the S = Q T = 1 pair correlation in the Bardeen-Cooper-Schrieffer (BCS) phase has been known 
to be very important [sol. |40|. From the analysis of the spatial structure of the Cooper pair in neutron matter, a size 
shrinking is found to occur at finite low density (2^ . [27l . [28j . It indicates that the spatially correlated neutron pair may 
appear in the transition region from weak coupling BCS regime to strong coupling BEC regime. Similar BCS-BEC 
crossover features of the spatially correlated neutron pair were predicted also at the surface of medium-heavy nuclei 
[13, mi- Since nuclear pairing increases as the size of nuclei shrinks (A ~ 12 /\fA MeV), one may expect that the 
BEC features are enhanced with respect to the infinite matter case. 

Thus, the spatial structure of neutron pairs at the surface of neutron-rich nuclei attracts presently great interest. 
In particular, BCS-BEC cross over features are expected to appear in the environment of dilute neutron matter, such 
as realized in neutron skins [U [12, [il| • This has been put into evidence in recent studies [13, [2l| where it was found 
that the extension of a nn-Cooper pair strongly decreases when going from inside to the surface before expanding 
again when leaving the nucleus. The minimum is very pronounced and the Cooper pair attains a minimum size of 
about 2-3fm. This is very small size reaching the dimension of the deuteron, that is of a bound state. This is also 
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highlighted by the fact that in -'^ Li even the single Cooper pair is very small in the surface [44|. 

The fact that ?in-Cooper pairs resemble a bound state in the nuclear surface and even become appreciably smaller 
than in low density neutron matter with maximum gap value, is intriguing. In this work we want to explore the 
rather extreme hypothesis that the binding mechanism may be influenced by some quasi two-dimensional features. 
By 'quasi two-dimensional (2D)', we mean that the 'bound' Cooper pairs are confined within a surface layer of about 
1 — 2 fm and that the degree of freedom in radial direction is approximately frozen. Such a scenario could for example 
be realized by the fact that the density distribution of valence neutrons in very neutron rich and heavy nuclei are 
in radial, say z-direction, concentrated in a surface layer and that the density distribution in z-direction can be 
approximated by a frozen Gaussian packet, whereas the motion of the neutrons within the layer is free. Pursuing the 
picture to its extreme, one could imagine a slab of low density neutrons where in the transverse (z) direction only 
a single Os level below the Fermi energy is active. Such situation also could eventually be relevant in the Lasagne 
phase in neutron star crusts, where a low-density neutron gas may appear as a layer on the surface of the Lasagne 
structure [i^ [H, [i^, SI] . In any case, one could imagine that in the outer layer of a neutron-rich nucleus, the pairing 
properties of the valence neutrons are more or less decoupled from the inert core part. For example, very neutron-rich 
Sn isotopes with > 82 could represent such a scenario. A 'quasi-2D' slab of low density neutron matter with 
the transverse degrees of freedom approximately frozen could then, as already mentioned, mimic the situation. We 
here insist on the 2D character because it is well known that the dimensionality has great influence on bound state 
formation and, thus, also on the formation of tight Cooper pairs. Such quasi-2D paired systems with low density 
also have been investigated in cold atom physics e.g. hydrogen atoms on a liquid "^He surface and laser cooled atoms 
trapped in an planar potential 0, [s^l • 

At any rate, and for the different reasons mentioned above, it seemed to us interesting to study an idealized situation 
of quasi-2D low density neutron matter realized within an infinitely extended slab. The study may turn out academic 
in the end but if future more detailed and not completely trivial investigations will exhibit a link to the surprisingly 
strong bound state features of Cooper pairs in the nuclear surface, mentioned above, our effort will have been worth 
while. In this respect it is of interest to mention that, indeed, our study will show that in 2D matter a Cooper pair 
attains a minimal size of about 2 fm whereas in 3D this value is over twice as large. 

In the present work, we investigate properties of neutron pairs in quasi-2D neutron systems by using finite-range 
effective nuclear forces. Concerning the degree of freedom in the transverse z direction perpendicular to the 2D plane, 
we for simplicity assume a Gaussian packet as mentioned above. We study the two-neutron wave function in this 
model of quasi-2D system and show that two neutrons form a tightly bound state with a small size when the width 
a of the Gaussian in z-direction is small enough. We also investigate pairing properties of neutron matter based on 
the BCS theory. Pairing gap and size of the Cooper pair in quasi-2D neutron systems are analyzed, and BCS-BEC 
crossover phenomena are discussed. 

The behaviour of T = 1 pairing in infinite three-dimensional (3D) matter has been extensively investigated based 
on the BCS theory [51|] for symmetric matter and neutron matter (for example, Refs. [13, [H, [2^ and references 
therein). As for a slab of nuclear matter, properties of T = 1 pairing in symmetric nuclear matter have been studied 
in Refs. [13, HI, [H, [IFl and it was shown that the pairing gap increases at the surface of the slab. As for T = 
pairing, the gradual transition from BCS to BEG involving the deuteron formation was predicted in dilute symmetric 
matter [13, [13, [13, [IE [Hsj • Competition of alpha and deuteron condensation was also discussed [13| ■ 

There are many varieties of BCS calculations of nuclear matter concerning, e.g., the treatment of single-particle 
spectra and the pairing interaction. For convenience, effective nuclear forces have been used in some calculations on 
pairing in nuclear matter instead of dealing with realistic bare nuclear forces. It was shown that the Gogny force 
|57ll58|. a finite-range effective nuclear force, gives similar behaviour of nuclear matter pairing to that calculated from 
realistic nuclear forces with the Brueckner-HF (BHF) approximation, as well as reasonable pairing properties in finite 
nuclei (2^. [stI. [ssI. [sqI. [60t . Therefore, we adopt the Gogny DIS force in the present BCS calculations. We also use 
the Minnesota force [61| which reproduces features of S-wave N-N scattering in T = and T = 1 channels, and 
nucleus-nucleus scattering in the light nuclear region such as the a-a system, as well as properties of light nuclei. The 
Minnesota force is often used in structure studies of light neutron-rich nuclei. 

This paper is organized as follows. In the next section, we explain our simplified model of quasi-2D systems and 
formulation of the present work. The effective nuclear forces are described in section Hill In section ITVl the results 
obtained for a two- neutron system and neutron matter are shown, and in section fVl discussions are outlined. Finally, 
we give a summary and outlook in section IVTl 

II. FORMULATION OF NUCLEAR 2D PAIRING 

Here we explain our framework of the simplified model of quasi-2D neutron system confined in a slab with a certain 
thickness. We first describe the wave function and the Hamiltonian for a quasi-2D neutron matter system. We then 
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explain the isolated two-neutron system and give the formalism for the quasi-2D system within the BCS theory. 



A. Quasi-2D neutron systems 



According to the outline of the Introduction, we propose a model for a system of neutrons confined in a slab with 
a certain thickness. We assume that the neutrons in transverse (z) direction are represented by a simple Gaussian 
packet of size a. As already mentioned, this size may mimic the concentration of the amplitudes of single particle 
wave functions in a surface layer, where z corresponds to the radial direction in finite neutron-rich nuclei. 

Therefore, the wave function for a single neutron is written as 



<i>(r) = $ 
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where indicates the coordinates in the slab, {x,y), and x is the intrinsic spin. The thickness, i.e., the rms width, 
of the slab is 2-\/ (z^) = \/2a. In a similar way, a iV- neutron wave function is given as 



'^(ri. 



i-N,Xi, ■■■,Xn) = *^^(ri±, • • • , tatj^, Xi, • • • , Xiv) ® <l>'"{zi) ■ ■ ■ <t>"'{zN)- 



(3) 



Here $ and are the antisymmetrized many-body wave functions. 

Neutrons are interacting with each other via two-body nuclear forces. The interaction V{rij) between two neutrons 
is the two-body effective nuclear force in the T — \ channel. Here r^j jri — rj|. Because the S'-wave neutron pair 
correlation is studied in the present work, only the central part of the effective nuclear force is considered. The Gogny 
and Minnesota forces are represented by a superposition of two Gaussians, 



V{r) = + - - MmPaPr) CXp 
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We freeze the the degree of freedom in transverse (z) direction, as described in the Introduction. By integrating 
over Zi coordinates the three-dimensional Schrodinger equation for $ is reduced to the two-dimensional equation for 
<i>^^ with respect to r^j^. Then we get the following equation for quasi-2D neutron systems, 



■ ,r7V_L-Xi, ■ 



,Xn) 



where the 2D-kinetic term and the 2D-interaction term are 



t± 



2mdY\ 2m \dx^ dy 



>rAr±,Xl7 



,XArj 



(5) 
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(7) 



y2^(r^) = {4P'{z,)<tP^{z^)\V{r)\4P'{z,)cp^^{z2)) 



J2 {a) {W„, + B,„P, - HmPr - M„,P„Pr ) CXp 
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Zb„ (a) 



62 + 2a2 
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(8) 
(9) 



As seen from ([8]) and ([9]), the strength of the two-body interaction V^^{r±) in the quasi-2D system depends on 
the thickness of the slab with the factor Zb^{a). The factor decreases with the increase of the width parameter a. 
The thickness of the slab is only reflected in a weakening of the interaction for bigger slab sizes in our genuinely 
two-dimensional model. These are, of course, strongly simplified, may be oversimplified, assumptions to describe the 
detailed behaviour of the valence neutrons because the Pauli effect from the core should be more complicated and a 
long tail of weakly bound neutrons may be not so simple. However, it is expected that such effects may effectively be 
taken into account by a modification of the width parameter a. In the present work, we analyze properties of neutron 
pairs in quasi-2D neutron systems for various values of the slab size a. 
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B. The isolated two- neutron system 



For later comparison we consider the S'-wave relative wave function Xi, X2) for an isolated two-neutron 

system in the quasi-2D space. It is obtained by solving the Schrodinger equation, 



I I d_ 

m \dr± r± dr± 



(10) 
(11) 



V'""(r_L) is the spatial wave function normalized as / (Pr±\iJ;™{rj_)\'^ = 1, and ^5=0 is the intrinsic-spin wave 
function for the spin-zero channel. 

It is important to note that in two dimensions, a bound state can form at any arbitrarily small attraction [63 |. 
Since the nuclear force in the channel is attractive at low energy, a two- neutron bound state in the quasi-2D space 
is obtained with an arbitrary width a. We calculate the exact solution of the bound state with a Gaussian expansion 
method. The size of the two-neutron bound state is given as 



d2rj.r^|?A""(ri)|2 



(12) 



Please note that the size in 3D is given as ^ = V (C"")^ + o,^- 



C. Quasi two-dimensional neutron matter within BCS theory 



In this section, we investigate the behaviour of neutron pairs in quasi-2D infinite neutron matter in the slab by means 
of the BCS theory, which is equivalent to the HFB approximation in a homogeneous case [s^llO]. The application of 
BCS theory with the Gogny force to a slab of nuclear matter was, e.g., already performed in Ref. [53| . 

Let us explain the equations for pairing in quasi-2D neutron matter. We apply the BCS theory to the 2D 
Schrodinger equation ((S]). The S'-wave pairing gap depends on |pj^| and it is written as 



A{p) 



1 f d^k. 



(2^) 



2«Pp(P-L - k_L) 



A{k) 



V(e(fc)-A^)2 + A(A)2' 



(13) 



Here p and k denote |p^| and |kx|, respectively, and Vpp is the pairing force, 

2 

^^Pp(P-L - ki) = ^ Zb^ (a) {W„i - B„i - Hrn + Mm) {\^brn) exp 



(14) 



(15) 



V (p) is the Hartree-Fock potential, 



V^'^ip) = f I] ZbM (2W^ + S„ - 2i/„ - M^) {V^bmf - / j^vp^ip^ - k 

m=l ^ ^' 



where the exchange term is given as 

2 



^'ph(Pi - k^) = ^ Zb,^{a) (Wm + 2S„ - H,n - ^M^n) {v^brn) CXp 



771 — 1 



^-bUp^-k^f 



The occupation probability is 



e(fc) — fi 



V(e(fc)-M)' + A(fc)2 



(16) 



(17) 



(18) 
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and the chemical potential /i is determined so as to satisfy the number equation, 

'4 f d-k^^^,^ ^^g^ 



p kp f (Pk± 
2 " 4^ Jq^"'" 



where p and kp are the density and the Fermi momentum in 2D, respectively. In this paper, the kp dependence of 
pairing properties in neutron matter is discussed as a function of the ratio kp/ko, where fco = 1.36 fm"-'^ is the Fermi 
momentum at the normal density of 3D symmetric nuclear matter. 
The energy per nucleon for the BCS state is 



Ebcs = -{'^>BCS\H^^\^BCS) = I „ 

p pJ [^T^r 

where the integrated pair energy Epair is 



2Wfe + EJ,a^r, (20) 



with 



^ — / ^A(fc)K(fc), (21) 
pJ (27r)^ 



1 A(fc) , , 

nik) — UkVk — — — , ^=. (22) 

^' 2 V(e(fc)-M)' + A(fc)2 ^ ' 



The energy per nucleon for the HF state $hf is determined from 

£;Hi. = i($HF|i?'^|$HF). (23) 
P 

To analyze the spatial structure of a nn-Cooper pair it is useful to study the pair wave function in the coordinate 
space. It is defined by the Fourier transform of the anomalous density K.{k), 

vl/p,i,(r^) ^ no($Bcs|anrY t)«^(rT i)|$BCs) ^ j ^n{k)e^'^-<^'--^'^\ (24) 

Here r±_ = \y'^_ — r" |, and no is the normalization factor so that J d^rx |5'pair(''±)P = 1- The pair size (coherence 
length) is calculated from the root-mean-square distance of the pair wave function 



where 

K)^ I d\^rl\^^.^,{r^)\\ (26) 



In the low-density limit kf ^ 0, the pair wave function ^'pair and the energy 2Ebcs approach -0"" and E'^^ of the 
two-neutron bound state of section II. B. 



III. EFFECTIVE NUCLEAR FORCES AND PAIRING PROPERTIES IN 3D NEUTRON MATTER 

As already mentioned, in this work we use the Gogny DIS force [l^jllll and the Minnesota force [6l|. In the T = 1 
channel, the Gogny force is density independent and is represented by a superposition of two Gaussians (Eq. 
The parameters for the Gogny DIS force are, 

bi = 0.7fm, 62 = 1.2fm 
Wi = -1720.30MeV, Bi = 1300.00MeV 
Hi = -1813.53MeV, Mi = 1397.60MeV 
1^2 = 103.64MeV, B2 = -163.48MeV 

H2 = 162.81MeV, M2 = -223.93MeV. (27) 
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The Minnesota force is also given by a sum of finite-range Gaussians, 

V{r) = ^{l-P„)^{l+uPr){VRCxp{^KRr^) + VsCxp{-Ky)) 

+ ^(1 + P^)\i'^ + ^Pr) {Vr exp{-KRr^) + Vt expi^Ktr^)) , (28) 

Vr = 200.0MeV,Kfl = 1.487fm"^ 
V, = -91.85MeV,K, = 0.465fm"^ 

Vt = -178.0MeV,Kt = 0.639fm"^ (29) 

We adopt for the exchange- mixture parameter the value of m = 1.0. Then the first (second) term gives the two-body 
force in the S* = T = 1 (S* = 1 T = 0) channel. 

In principle the Gogny force represents an effective in medium force and the Minnesota force a phenomenological 
nuclear force in free space. It has, however, been shown that the Gogny DIS force gives similar behaviour for 
the pairing gap of symmetric nuclear matter as that calculated from realistic nuclear forces with the BHF+BCS 
approximation, and also describes pairing properties in finite nuclei [U, [13, [H, [H, [13] • The Minnesota force reasonably 
reproduces the S'-wave N-N scattering in the S = T = 1 channel as well as the S = 1 T = channel. It also 
well describes nucleus-nucleus scattering such as the a-a system, as well as properties of light nuclei. Because the 
Minnesota force has no density-dependent term in the S = 1 T = channel, it can not describe the saturation 
properties of symmetric nuclear matter. However, it gives a similar density dependence of neutron matter energy as 
the one obtained with the Gogny DIS force, as shown later. 

The values of the n-n scattering length in the channel obtained with the Gogny DIS force and Minnesota force 
are Us ~ —12.1 fm and as = —16.8 fm, respectively. These values semiquantitatively agree with the experimental value 
fls = —18.5 ±0.4 fm 'so]. In Fig. [1] the calculated phase shifts of N-N scattering in the channel are shown as well 
as those obtained with a realistic force. Both the Gogny DIS and Minnesota forces show a reasonable reproduction 
of the ^5*- wave scattering. 

Figure shows the properties of infinite 3D neutron matter calculated with the Gogny DIS and Minnesota forces 
by means of the BCS theory as well as from the HF calculations. Figures UK a) andHJb) show the the energy {Ercs) 
and the integrated pair energy {Epair) per nucleon, and the chemical potential (/z) of the BCS state are shown as well 
as the energy of the HF state (Ehf) for 3D neutron matter with the Gogny DIS and Minnesota forces, respectively. 
The density dependence of the neutron matter energy obtained with both forces is in reasonable agreement with that 
calculated with ab initio methods using realistic nuclear forces [H, . 

In Fig. EJc), the pairing gap of neutron matter is shown. The difference between the Gogny DIS and Minnesota 
forces is not large in the low-density region kp/ko < OA. In the kp/ko > 0.4 region, the Minnesota force gives a much 
smaller pairing gap than the Gogny DIS force because it has a larger repulsive core with 108.15 MeV hight at r = 0. 

Let us make a comment on the pairing gap in 3D neutron matter investigated in different approaches. There are 
many calculations of the pairing gap in nuclear matter using realistic nuclear forces (for reviews [H, [Hg). In the 
standard BHF+BCS approximation, the pairing gap in 3D neutron matter is about 3 MeV at the peak position 
kf/ko ^ 0.6 (figure 12 of Ref. [1^), which is consistent with the pairing gap calculated in BCS theory with the Gogny 
force. In different approximations with medium polarization effects, a quenching of the pairing gap in neutron matter 
was suggested [M [Hi [H [13, [H, [H . 

The size ^ of a Cooper pair in 3D neutron matter is plotted as a function of kp/ko in Fig. [D(d). As is already 
discussed in Refs. [23. [27l [28|. the pair size shrinking in 3D neutron matter is seen at finite low kp. Similarly to the 
pairing gap, the difference of the pair size between the Gogny and Minnesota forces is small in the low-density region. 
In the kp/ko > 0.4 region, the pair size increases more rapidly because of the weaker pairing with the Minnesota force 
than is the case with the Gogny DIS force. As we show later, also in quasi-2D neutron matter these two forces give 
qualitatively similar pairing behaviour for low kp. With both forces, the ^ < d region [d is the average inter-neutron 
distance d = p~^^^) appears at finite low kp, i.e., at low density. It shows the BCS-BEC crossover feature as suggested 
by Matsuo (27| . 

IV. RESULTS 

In this section, two-neutron correlations in quasi-2D neutron systems are investigated by performing BCS calcu- 
lations for 2D neutron matter with the Gogny DIS and Minnesota forces. The results for various values of a of the 
confinement in the frozen direction (z) are analyzed, and the behaviour of the Cooper pair in quasi-2D neutron matter 
is discussed as a function of kp/kg. 
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FIG. 1: Phase shift of A^'-A'' scattering in the channel obtained with the Gogny DIS force and Minnesota force. phase 
shift of a realistic force, Reid93 potential [70|], which fits the experimental proton-proton scattering cross sections, is also shown. 



A. Two-neutron bound state in quasi-2D system 



By solving the Schrodinger equation pTI)l . the wave function (r^) of two isolated neutrons in a quasi-2D system 
is obtained. While a two-body bound state in 3D forms only at sufficiently strong attraction, a bound state in 2D 
can form, as already mentioned, at any arbitrarily small attraction. Since the nuclear force is attractive at low 
energy, the two-neutron bound state is formed in quasi-2D systems with arbitrary width a. 

The energy E^^ and the pair size ^""^ of the bound state for a = 1 , a = 2 and a = 4 fm are given in Table [H 
Energy and size of the dineutron decrease as the width a becomes small. In the quasi-2D system with a small width 
of a = 1 fm, the dineutron size is as small as that of the deuteron in 3D. On the other hand, in the quasi-2D system 
with a width a = 4 fm, the dineutron has a small binding energy, less than 0.1 MeV, and a size larger than 20 fm. 
We, therefore, see that in our model the width of the slab has a strong influence on the two body correlations. 

It is very interesting that, although two neutrons are unbound in 3D, they form a tightly bound dineutron state with 
a size comparable to deuteron dimensions when they are confined in a sufhciently thin slab. To form the small-sized 
dineutron, the reduction of the dimension from three to two is crucial. In the present model, neutrons are confined 
in a Gaussian packet with respect to the frozen ^-direction. 

In slab direction, the potential part V^^ {r± ) is given by almost the same form as the original effective nuclear force 
V{r) except for the factor Zb^{a) in ([8|). This factor is smaller than one, depending on a. For example, for the case 
of a = 1 fm and the interaction range 62 = 1.2 fm as with the long-range attractive term of the Gogny DIS force, 
one gets Zf,^(a) = 0.59, that is the strength of F^^(r_L) is smaller than that in the original interaction V(r). Let 
us remind the reader that the size of the dineutron state is determined by the balance of the kinetic energy and the 
potential energy. In a 2D system, the kinetic effect is weaker and a small-sized two-particle bound state can more 
easily be formed than in the 3D case. 

In Fig. [31 the wave function ?/;""(rj_) and the probability density r'_L|?/;""(rx)p of two isolated neutrons in a quasi- 
2D system are shown as a function of the relative coordinate r^. In the quasi-2D system with a = 1 and a = 2 fm, 
the probability density has a peak around ri_ — 2 fm, which shows the significant spatial correlation between the 
neutrons. In a quasi-2D system with a = A fm, the density distributes in a broad region, and the spatial correlation 
is much weaker. 



B. Quasi-2D neutron matter 



We now extend the BCS calculations to quasi-2D neutron matter with varying width parameters a. 

The pairing gap at the Fermi momentum A^;- = A(fc = kp) as a function of kp/ko is shown in Fig. [4l The Gogny 
DIS and Minnesota forces give qualitatively similar results in the low kp region. At kp/k^ > 0.3, the pairing gap 
is more suppressed in case of the Minnesota force, because, as already mentioned, the short-range repulsion is larger 
than with the Gogny force. For a = I fm, the pairing gap obtained with the Gogny DIS force has a peak about 4 
MeV high at kp/kg ~ 0.4, while that calculated with the Minnesota force has a peak ~ 3 MeV high at kp/ko ~ 0.3. 
With the increase of the width a of the slab, the pairing gap Aj? is quenched because the attraction of the pairing 
force is weaker in a wider slab. 
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FIG. 2: Energy, pairing gap, chemical potential and pair size in 3D neutron matter as a function of kp/ko {ko = 1.36 fm~^) 
obtained with the Gogny DIS force and Minnesota force. (a)(b) Energy (Ebcs), the integrated pair energy (Epair) per nucleon 
and the chemical potential {^) of the BCS state, and the energy of the HF state (Ehf)- (c) Gap energy at the Fermi momentum 
Af = A(fc — kp)- (d) Pair size ^. The average inter-neutron distance d = p^^^^ is shown by a thin dashed line. 



In Fig. [51 the energy {Ebcs) and the integrated pair energy {Epair) per nucleon, and the chemical potential (/i) of 
the BCS state are shown as well as the energy of the HF state {Ehf) calculated without pairing. Compared with the 
energy of 3D neutron matter shown in Fig. [21 one of the striking features of quasi-2D neutron matter is that EbcSi 
Epair, and the chemical potential have finite negative values even in fcp — > limit. This is because, in the kp 
limit, the system goes to the strong coupling BEC limit, where Cooper pairs become identical with the two-particle 
bound state. In this limit, —2Ebcs and —2fi become equal to the dineutron binding energy. As soon as the chemical 
potential turns negative, there is true binding between the two neutrons and the strong coupling BEC feature is born 
out. Therefore, in quasi-2D neutron matter one can, contrary to the 3D case, follow the transition from BEC to the 
BCS regime. On the other hand, because of the same reason, the pairing gap Ap in the quasi-2D neutron matter 
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TABLE I: Energy El° and size 5"" of the two- neutron bound state -(/^""(rx) in the quasi-2D system with the width parameter 
a = 1, 2, and 4 fm. The energy and the size for a deuteron in 3D space obtained with the Minnesota force is also shown for 
comparison. 



Gogny DIS 
a = l a = 2 a = 4 
El'^ (MeV) -2.7 -0.64 -0.03 

(fm) 3.5 6.8 33 

Minnesota 
a=l a = 2 a = 4 
El° (MeV) -1.9 -0.74 -0.06 
^i" (fm) 4.3 6.5 21 

Deuteron with Minnesota 
E (MeV) -2.2 
C^" (fm) 3^} 
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FIG. 3: Two-neutron wave function ^""(rx) in quasi-2D systems with the width parameter o = 1, 2, and 4 fm calculated with 
(a) the Gogny DIS force and (b) Minnesota force. (c)(d) Same as (a)(b) but for the probability density rxiV'""(''-L)P- 



shows a linear behaviour at the low kp limit (Fig. 2]), while the pairing gap disappears with a horizontal tangent at 
small in 3D neutron matter (Fig. [2l^c)). 

The pair size in quasi-2D neutron matter is plotted as a function of kp/ko in Fig. [6l In the kp limit, 
the pair size equals the size of the two-neutron bound state (^2" '^^ CH))- With the increase of kp, the pair size 
reduces first before it becomes large again with further increase of fc^^. Thus the minimum pair size is found to be 
smaller than the size of the isolated two-neutron bound state in the quasi-2D system. This size reduction of the 
nn-Cooper pair indicates the enhancement of the dineutron correlation at finite low kp, i.e. due to the existence of a 
Fermi sea. It should be pointed out that the pair size significantly depends also on the width a of the slab. 

A striking difference of the pair size between the quasi-2D and 3D neutron matter is that the pair size at fci? 
limit is finite in the quasi-2D neutron matter while it becomes infinite in the 3D neutron matter. In the 3D neutron 
matter, the Cooper pair size ^ has a minimum around kp/k^ — 0.5, and it rapidly increases as kp goes 0. Interestingly, 
the size shrinking of the Cooper pair in quasi-2D neutron matter has a close analogy with that of the deuteron in 3D 
symmetric matter predicted by Lombardo and Schuck [s^l . 

As mentioned above, the spatial correlation of Cooper pairs enhances at finite low kp in quasi-2D neutron matter 
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FIG. 4: Pairing gap at the Fermi momentum Af — A{k — kp) in quasi-2D neutron matter as a function of kp/ko {ko — 1.36 
fm^^). Upper and lower panels show the results with the Gogny DIS force and Minnesota force, respectively. 





FIG. 5: Energy (Ebcs), the integrated pair energy (Epair) per nucleon and the chemical potential (p) of the BCS state, and 
energy of the HF state (Ehf) as a function of kp/ko {ko = 1.36 fm~^). (a) and (b): Energy of the quasi-2D neutron matter 
with a = 1 fm obtained with the Gogny DIS force and Minnesota force, (c) and (d): Energy of the quasi-2D neutron matter 
with a = 2 fm. 



with a small width a, where the BCS-BEC crossover phenoniena is expected. To discuss pairing properties from 
the point of view of the BCS-BEC crossover, it is useful to compare the pair size (^±_ with the average inter-neutron 
distance d = p~^^^ as discussed in the works for 3D neutron matter [27,. i28j. The ratio (,±/d in quasi-2D neutron 
matter is shown in Fig. [7|as well as ^/d for the 3D neutron matter. Strictly speaking, because the transition from the 
weak coupling BCS phase to strong coupling BEC is known to be a crossover [t^, there is no critical boundary of the 
transition. However, we consider that BCS-BEC crossover features arise when the Cooper pair size is smaller than the 
average inter-neutron distance by adopting the same criterion as in Refs. [l^, H^. With this criterion, the BCS-BEC 
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FIG. 6: (a)(b) Size of a Cooper pair \l/pair in quasi-2D neutron matter obtained with the Gogny DIS and Minnesota forces. 
The calculated values are plotted as a function of kp/ko {ko = 1.36 fm~^). The average inter-neutron distance d = 
plotted with the dotted lines. (a')(b') Same as (a)(b) but scaled up for the vertical axis. 



crossover phenomena appears at kp/ko < 0.4 ^ 0.5 in quasi-2D neutron matter with a = 1 fm, and kp/ko < 0.3 ~ 0.4 
for a = 2 fm, while in 3D neutron matter it is seen at kp/ko < 0.4 ^ 0.5. In the quasi-2D neutron matter, the ratio 
£,±/d monotonically deceases a.s kp becomes small, and finally the system goes to the strong coupling BEC limit. 

To reveal the features of BCS-BEC crossover it is also useful to analyze the spatial structure of the nn-Cooper pair 
wave function [89| as done in Refs. [13, HI]. Figures [S] and [H] show the probability density r±\^paiT{r±)\^ and pair wave 
function ^pa,ir{r±) of the Cooper pair in the quasi-2D neutron matter. Let us discuss the results for a = 2 fm. At low 
kp such as kp — ko/8 and kp = fco/4 (Fig. [51^b4) and (b3)), the spatial correlation of the Cooper pair is strong and it 
indicates a BEC-like dineutron pair. From kp — ko/8 to kp = /eo/4, the first peak becomes narrow due to the Pauli 
principle, stemming from the other Cooper pairs. Consequently, the size shrinking of the dineutron occurs at a finite 
low kp, as mentioned. With further increase of kp, the nodal structure appears at kp = fco/2 which indicates the 
transition from the BCS-BEC crossover to the BCS regime, and at kp = ko the short-range correlation disappears as 
in the weak coupling BCS phase. 

The pair wave functions in quasi-2D neutron matter with a = 2 fm are shown in Fig. [5) We also show the two- 
neutron wave function il]™(r±/) for the dineutron bound state in the quasi-2D system and 5'pair(7') in 3D neutron 
matter. In the kp ^ Q limit, the pair wave function ^'pair(''j_) becomes equal to the isolated two-ncutron wave 
function tjj™{r±) of the dineutron state in case of quasi-2D neutron matter, while ^'paii ('') in 3D neutron matter 
becomes infinitely broad because there is no two-neutron bound state in a 3D system. 



V. DISCUSSION 



A. Possible relation to neutron pair correlations in neutron skins 



In this section we tentatively associate the strong pairing correlations of valence neutrons in neutron-rich nuclei 
with the ones in a 2D system. It has recently become evident that in neutron-rich or neutron halo nuclei pairing 
correlations may show a strong maximum in the surface. This is for instance evidenced by the fact that the size of 
the Cooper pairs has a sharp minimum in the surface. This feature has been shown to hold true in the case of a single 
Cooper pair [tJ as well as in cases of many Cooper pairs f^H. Though for the moment, it cannot be excluded that 
the pronounced small Cooper pair sizes in the nuclear surface is due to a somewhat trivial, i.e. geometrical effect fTSj. 
we here make the hypothesis that this is due to a real strong correlation effect, due to the presence of a surface and 
hence of a 2D influence. 

The valence nucleons may be confined in a surface shell for different reasons. For example, if valence neutrons 
are, in a heavy nucleus, distributed around a high lying principal shell, the corresponding density distribution, for 
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FIG. 7: Ratio of the nn-Cooper pair size to the average inter-neutron distance d as a function of kp/ko {ko — 1.36 fm~^). 
(a)(b) Ratio ^±/d [d — p~^^'^) in quasi-2D neutron matter obtained with the Gogny DIS and Minnesota forces. The /cf^x ~ 1 
hne for a typical value of the BEC-BCS crossover [tH. [t^. ItsI] in 2D is shown by solid thin lines, (c) Ratio ^/d {d — p^^^^) in 
3D neutron matter with the Gogny DIS and Minnesota forces. The short-dashed lines indicate the boundaries for the weak 
coupling BCS region $ , /d > 1.1, the crossover region 0.19 < ^/d < 1.1 and the strong coupling BEC region ^/d < 0.19 in the 
regularized model [2?!. ItI. [75(| . Thin solid line is ^/d = 0.36 for the unitary limit. 



instance if the shell does not contain an S wave, may be concentrated in the surface region. Even though low angular 
momentum orbits with radial nodes have some strength distribution in the inner region, their probability amplitudes 
may show peak structure in the surface, and therefore, the pairing is expected to be induced by the surface neutrons 
due to the possible mixing of many states having significant overlap of their radial part in the surface region. 

On the other hand it is also well known that finite size enhances pairing. This is already documented by the fact 
that the nuclear gap varies on average like 12/\/A MeV, that is the presence of a surface certainly plays a role. This 
finite size enhancement is a quite subtle effect as theoretical investigations show (79j and implies several competing 
effects. How much the scenario is influenced by quasi 2D features is still an open question but an interesting and 
suggestive hypothesis which we wanted to explore in this work. So one may consider that the strongly correlated 
Cooper pairs with about a rms radius of 2 fm live in a surface sheet about 1 — 2 fm across. In section |ll] we mimicked 
the transverse motion of the neutrons by a simple frozen Gaussian (jjoi^). The motion perpendicular to z is considered 
as planar. This is, of course, a highly idealized representation of the situation as mentioned several times. It may 
nevertheless capture some 2D features which possibly underly the strong pairing correlations in the neutron skin of 
neutron-rich nuclei. 

In the previous analysis of pair size in quasi 2D neutron matter in section HVl the BCS-BEC crossover may occur 
only in dilute 2D matter for densities p < 0.05 — 0.07 fm^^. Let us estimate how many small sized Cooper pairs can 
condensate in an area AwR^ of the nuclear surface. For the dineutron condensation, more than one pair should exist 
at least. From the condition p < 0.05 — 0.07 fm^^, two dineutrons can condensate in the surface area for i? > 2 — 3 
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FIG. 8: Probability density r^|^'pair(»'x)P of a nn-Cooper pair in the quasi-2D neutron matter as a function of the relative 
distance r±. (al)(a2)(a3)(a4) Probability density in the slab with a = 1 fm at fc^ = fco, kp ~ ko/2, kp = fco/4, kp = fco/8. 
(bl)(b2)(b3)(b4) The density in the slab with a = 2 fm. (cl)(c2)(c3)(c4) The density in the slab with a = 4 fm. 



fm, and four dineutrons for i? > 3 — 4 fm. 

Recently, two-neutron correlations in neutron-skin nuclei have been investigated in light nuclei such as ^He and 
^^C [l3, [11, as well as medium-heavy nuclei (13, HH- *He is considered to be a neutron-skin nucleus having four 
valence neutrons surrounding an a core, and therefore this nucleus is one of the simple examples to discuss neutron 
pair correlations in neutron skins. In Refs. [13, [H, [3l, it was argued that the ^He ground state has spin-zero two- 
neutron correlations. From the above-mentioned discussion of quasi-2D neutron matter, the correlation in the ground 
state is expected to be of BCS-type rather than in the BCS-BEC crossover region, because the valence neutron 
distribution concentrates in the surface with a distance R < 2 fm, which is smaller than the BCS-BEC crossover 
condition i? > 2 — 3 fm for the condensation of two dineutrons. Instead of the ground state, the possible existence of 
an excited state ^He(Oj) with a dilute gas-hke state of dineutrons was suggested by one of the authors [l^l- According 
to the theoretical results of Rcf. 17], the dineutrons exist around the region with i? = 4 — 5 fm which satisfies the 
above-mentioned condition for the BCS-BEC crossover. Therefore, this state could have BEC-like dineutrons at the 
surface. 

In medium-heavy nuclei, the spatial structure of spin-zero neutron pairing at the surface was investigated based 
on HFB calculations in Refs. [20, 21]. To discuss two-neutron correlations in the neutron skin in connection with 
quasi-2D neutron matter, one should analyze the radial and angular behaviour separately and focus on the angular 
correlations (opening angle) of neutron pairing, as well as on the size of the Cooper pairs. The transition from BCS 
phase to BCS-BEC crossover is characterized by the change of the pair wave function from a weak-coupling nodal 
structure to an enhanced single bound state like peak. If such a structure change is found in the angular behaviour 
it may be an indication of the BCS-BEC crossover, while a clear evidence may not be seen in the radial structure 
because it might be affected by Pauli effect from constituent neutrons of the core. The spatial structure of spin-zero 
neutron pairing is displayed in Ref. [20| . which shows that the nodal structure along the angular direction is seen 
inside from the surface and it disappears on the outer side of the surface. More detailed analysis of the angular 
behaviour is requested. On the other hand, a good measure of the correlations also is the size of the Cooper pairs. 
Their rms radius seems systematically very small in the outer part of the nuclear surface [2lj, that is much smaller 
than the smallest value in homogeneous 3D matter. Whether this also always goes along with a maximum in the pair 
correlations at the same place in a finite nucleus, remains to be seen. 

Let us also comment on other effects omitted in the present model of infinite quasi-2D neutron matter. At the 
surface of a realistic finite nuclear system, valence neutrons feel the centrifugal force and the spin-orbit force which 
may work so as to break nn-Cooper pairs. Periodicity condition on the spherical surface is another effect to be 
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FIG. 9: (a)(b)(c)(d) Pair wave function \!'pair of a nn-Cooper pair in the quasi-2D neutron matter with a = 2 fm as a function 
of r±, and (f)(g)(h)(i) that in 3D neutron matter as a function of r. (e) Two-neutron wave function il)""{r±) of the bound 
state in the quasi-2D system with a = 2 fm. 



considered. 



B. Relation to two- neutron halo nuclei 



Two-neutron correlations have extensively been investigated in two-neutron halo nuclei by means of three-body 
models of two neutrons and a core 0, i, 1 1 1 0, i, i 0, Cl 113, iH, HI . Recently, it was suggested that the internal 
two-neutron wave function analyzed as a function of the distance R from the nuclear center shows behaviour similar 
to the pair wave function of 3D neutron matter in the BCS-BEC crossover region fF^. In the analysis of Ref. ^1], the 
inter-neutron wave function has a nodal structure at small i?, and it shows enhancement of the spatial correlations 
with increasing R. Moreover, in a further analysis of the neutron pair in ^"'^Li, size shrinking was found at the surface 
p3 | . It was found that there the pair is as small as 3 fm at i? ~ 3 fm 'ii'l . The formation of the small sized dineutron 
at the surface of ^^Li may be understood with the picture of a quasi-2D neutron system. As shown in Ref. |10j] with 
a three-body model, the spin-zero component of two neutrons shows a strong angular correlation and its distribution 
concentrates at the surface with i? = 3 ~ 4 fm. The situation corresponds to the quasi-2D neutron where the valence 
neutrons are confined at the surface region and the neutron pair is formed due to the angular correlation. It may 
be interpreted as the reduction of the dimension from three to two, being crucial for the small-sized neutron pair in 
two-neutron halo nuclei. At the small R region, i.e. near the core, the size of the neutron pair increases and nodal 
structure appears in the relative wave function in the three-body model. This is because the Pauli effect from the 
constituent neutrons of ^Li core raises the Fermi momentum of the valence neutrons. This behaviour is similar to the 
weakening of the pair correlation in the BCS region at kp > ko/2 in the quasi-2D neutron matter. 



15 



In spite of the analogous behaviour of the neutron pair in ^^Li to that of the Cooper pair in the quasi-2D neutron 
matter, one should keep in mind that two-neutron halo nuclei has only one pair of neutrons and the phenomena there 
is not regarded as the Cooper pair condensation in the BCS picture. Several investigations have, however, revealed 
that the features of a single Cooper pair and the ones of a Cooper pair embedded in several others are not drastically 
different. 



C. Two-proton correlation in quasi-2D systems 

It is also interesting to consider two-proton correlation in proton-rich nuclei in analogy to the two-neutron correlation 
in neutron-rich nuclei. While the energy of infinite proton matter diverges without charge neutrality, proton-rich 
matter such as proton halo or proton skin structure can exist in finite nuclei. Moreover, since two-proton decays are 
allowed from some proton-rich nuclei such as ^Be,®C,^^0 and ^^Ne [H, [13, HE]) more direct information of T = 1 pair 
correlation could be experimentally obtained than for the case of neutron-rich nuclei. 

In the same way as neutron systems, we calculate a wave function 'ip'^"{r±) for the two-proton bound state in 
quasi-2D systems. Coulomb force (Vcoui(f)) is added to the two-body effective nuclear forces. In Fig.lTOl the obtained 
two- proton wave functions in the ^5- wave channel are shown in comparison with those of the dineutron state. Because 
of the Coulomb force, the spatial correlation of the T = 1 pair is slightly more suppressed in the two-proton system 
than in the two-neutron system. However, the wave functions for a = 1 fm and a = 2 fm still have a pronounced 
peak. It suggests that the spatial correlation of the T = 1 pair also is present in proton systems. 
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FIG. 10: Two-proton wave functions ip^''{r±) in quasi-2D systems with the width parameter a = 1 and a = 2 fm calculated 
with (a) the Gogny force and (b) Minnesota force complemented with Coulomb force (Vcoui). Two-neutron wave functions 
'4^""{r±) (without Vcoui) are also shown for comparison. 



VI. SUMMARY AND OUTLOOK 

Singlet S (^S) nn correlations in quasi-2D neutron systems were studied in a simplified model of neutrons confined 
in a slab with a certain thickness. The finite-range effective nuclear forces, the Gogny DIS and Minnesota forces, were 
adopted for the study. The strength of the two-body interaction in the quasi-2D systems depends on the width a of 
the slab. Namely, the interaction is stronger for thinner slabs than for wider ones. When the width of a slab is small 
enough, e.g. around a — 1 fm, two neutrons form a tightly bound dineutron with a size as small as the deuteron size. 
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The behaviour of the nn-Cooper pair in the quasi-2D neutron matter was investigated by means of the BCS theory. 
In quasi-2D neutron matter, the transition region of the BCS-BEC crossover appears at kp/ko < 0.4 ^ 0.5 and for 
a ^ 1 fm. The size shrinking of the Cooper pair occurs at finite low kp because of the PauU effect. It is also important 
to note that, as the width a of the slab decreases, the pair correlation in 2D is enhanced and the Cooper pair size 
becomes small because of the stronger attraction in the thinner slab. This is one of the characteristics different from 
3D. In particular, the reduction of the dimension from 3D to 2D is crucial to form the strongly correlated neutron 
pair. The chemical potential becomes negative for kp/ko < 0.2 and, therefore, there is true two-neutron binding in 
quasi-2D neutron matter and a genuine BEC regime is born out. 

The relation of the quasi-2D neutron matter to finite neutron-rich nuclei was also discussed. Several studies of 
the past have revealed that the size of Cooper pairs changes very substantially from inside to outside of a nucleus 
[21I [4^ . The subject has been taken up again more recently and the radius dependence of a neutron-neutron (nn) 
Cooper pair has been investigated more systematically [2l|, Q] . It thereby turned out that the old believe that the 
size of a Cooper pair is of the order of the diameter of a large nucleus must be revisited when considered locally, i.e., 
as a function of the nuclear radius. Namely, detailed BCS calculations with the DIS Gogny force have shown that 
the rms radius of an nn-Cooper pair, though indeed large in the interior, becomes very small in the outer surface of 
a nucleus, where the density is low, before expanding again when leaving the nucleus. The minimum extension of a 
Cooper pair attains a size as small as 2 fm, that is the same as the size of a deuteron, a bound state. Although the 
variation of the pair radius as a function of the density qualitatively follows the one of infinite matter, a factor of two 
reduction of the pair diameter with respect to the homogeneous case is a quite surprisingly strong effect which calls 
for an explanation. 

We considered this pair size shrinking in the surface of neutron-rich nuclei with the present calculations of quasi-2D 
neutron systems, based on the hypothesis that the valence neutrons in a neutron-rich nucleus live in radial direction in 
a surface sheet about 1 — 2 f m across because the wave functions of the active shell are essentially concentrated in the 
surface while the nucleons are free to move within the sheet. The idea behind this scenario is that the neutrons in the 
surface sheet are essentially constrained to a two-dimensional world which could give a certain clue to the smallness of 
the Cooper pairs, since in two dimensions even an isolated nn pair forms a bound state. Pursuing this idea, we further 
idealized the situation assuming a planar slab of neutrons where for simplicity the motion perpendicular to the slab 
is frozen into a Gaussian wave packet in order to simulate the surface concentration of the valence neutrons in radial 
direction. Freezing the perpendicular motion means that the neutrons are two-dimensional in the direction of the slab. 
The thickness of the slab only influences the strength of the pairing force, as already mentioned above. We, indeed, 
find that for reasonable slab thicknesses of a ~ 1 fm, the pair radius is of the order of 2—3 fm in agreement with 
realistic 3D HFB calculations [IJ. We are, of course, aware of the fact that our model may be strongly oversimplified 
and further studies are required in order to ascertain whether a 2D aspect indeed helps to interpret the occurrence of 
astoundingly small Cooper pairs in the nuclear surface. One also should keep in mind that a local picture may be a 
little too simple in the sense that the pair size certainly fluctuates a lot. Nevertheless, a classical picture often grasps 
a good deal of reality, even in so tiny quantum object as nuclei. Neutron pairs much smaller than the dimension of 
a nucleus may approximately be considered as bosons and a strong coupling Bose-Einstein condensation scenario of 
nn pairs may be realized in the surface of neutron-rich nuclei. Bound state like Cooper pairs in the nuclear surface 
may also allow to reinterpret neutron transfer reactions into superfluid nuclei which are known to have a strongly 
enhanced cross section j86j. 

2D aspects may also be relevant in the Lasagne phase in neutron star crust. In practical application to such systems, 
one should extend the model by taking into account finite temperature and also excitation in the z-direction for higher 
neutron density. 

BCS and BEC problems of two-dimensional fermion systems have been discussed also in the field of condensed 
matter physics in relation to electron systems and atomic systems [tiI, [T^I- Phenomena of BEC and BCS-BEC 
crossover in quasi-2D systems of trapped gases have been discussed in recent works [13, [13: [11|: and it has been 
suggested that the condensation in a quasi-2D gas is sensitive to the frequency u of the confinement in the "frozen" 
direction. This feature may have an analogy to the present results for the quasi-2D neutron matter. 
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